Critical state in type-II superconductors of arbitrary shape 
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The well-known Bean critical state equations in general are not sufficient to describe the critical 
state of type-Il superconductors when the sample shape is not symmetric. We show how one 
can find the critical state in superconductors of arbitrary shape. Analyzing a simple example of 
nonsymmetry, we demonstrate that in the general case, a perturbation of the current distribution 
in the critical state propagates into the sample smoothly in a diffusive way. This is in contrast to 
the usual Bean critical state where the current distribution changes abruptly at a narrow front. 
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The concept of the critical state introduced by Charles 
Beari is widely used to describe various physical phe- 
nomena in the vortex phase of type-II superconductors, 
see, e.g., Refs. 1m 

and citations therein. According to 
Bean, in the critical state of type-II superconductors with 
flux-line pinning, the driving force of the currents flowing 
in this state is balanced by the pinning force acting on 
the vortices. The critical state is characterized by the 
component of the current density flowing perpendicular 
to the flux lines, jc±) since only this component gener- 
ates a driving force. It is assumed in the critical-state 
theory that this Jc_l is known, i.e., it is a given func- 
tion of the magnetic induction B, jci_ — jc±(B), and the 
problem of this theory is to find the appropriate distri- 
bution of the magnetic fields and currents in the critical 
state. Below, to explain the physics with the least mathe- 
matical complications, we shall imply the simplest form: 
jc_L = const., which is frequently used in practice, but 
an extension to the general case is straightforward. For 
simplicity, we also assume that the magnetic fields H in 
the superconductor considerably exceed the lower critical 
field Hci, and so we may put B — /xqH. 

If in the critical state a current-density component j\\ 
parallel to the local magnetic field is also generated, the 
magnitude of j|| remains undefined, and one thus can- 
not in general find the distributions of the magnetic field 
H(r) and current density j(r) in the critical state. In- 
deed, to solve the Maxwell equations for H, 



rotH=j, divH = 0, 



(1) 



it is necessary to know the magnitude and direction of 
the currents j(r) in the sample.^ However, one has only 
the two conditions: 



j±=jc±, divj = 0, 



(2) 



for three quantities: j|| and the angle defining the 
direction of j± in the plane normal to H(r). Thus, the 
existing critical-state theory based only on Eqs. (0) 
is not complete. 

The above Maxwell equations with conditions ||2Jl can 
provide the description of the critical state when the 



shape of the superconductor is sufficiently symmetric and 
the external magnetic field is applied along a symmetry 
axis, so that some constraint on the directions of the cur- 
rents is known in advance. For example, the direction of 
the currents is obvious for a slab in an external magnetic 
field parallel to its surface. For an infinitely long cylinder 
with arbitrary cross-section in a magnetic field parallel to 
its axis, the currents flow perpendicular to this axis, and 
the critical state problem is solvedi^ Another completely 
solvable case is infinitely thin flat superconductors^ for 
which the currents can flow only in the plane. However, 
in the case, e.g., of a thin rectangular platelet oi finite 
thickness in a perpendicular magnetic field, the above 
critical state equations are already incomplete for deter- 
mining the magnetic fields and currents in the critical 
state.^ 

We emphasize that even for simple experimental situ- 
ations equations together with conditions l(2Jl can be 
insufficient for solving the critical-state problem. As an 
example that we shall analyze below, consider an infi- 
nite slab of thickness d. Let this slab fill the space |a;|, 
y| < oo, \z\ < d/2, and be in a constant and uniform 
external magnetic field Ha {Ha ^ Jc = jc±d) directed 
along the z axis, i.e., perpendicularly to the slab plane. 
Let then a constant field hax {Jc/'2 < hax ^ Ha) be ap- 
plied along the x axis, and after that the magnetic field 
hay {hay <Si Ha) Is switchcd on in the y direction. This 
critical state problem is not fully defined. Indeed, the 
condition divj = yields jz = 0, i.e., the currents flow in 
the x-y planes. Then, to describe the critical state, we 
may use the parametrization: 

j = ic{^,e,%l}){cosLp{z), sin(^(z),0), 
H(z) = H, + h(z), 
h(2) = {hx{z),hy{z),Q), 

where jc(<P, Q, tp) is the magnitude of the critical current 
density when a flux-line element is given by the angles ip 
and 61, tauT/i = hy/h^, tan9 = {h'l + hy)^/'^ / Ha, while the 
current flows in the direction defined by the angle tp; all 
these angles generally depend on z. A dependence of jc 
on the orientation of the local H, jc{f,0,ip) — ic±/[l ~ 
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cos^((^ — -0) sin^ 0]^^/^,^ appears if jc is not perpendicular 
to this H. However, at Ha ^ hax, hay, Jc, the field H is 
practically normal to the x-y planes where the currents 
flow {9 ~ 0), and we may put jc{f,0,^) = jc±- With 
this parametrization, the equation divH = is satisfied 
identically, while the Maxwell equation rotH = j reads 



the total current density j(r). Using the Maxwell equa- 
tion 



dz 
dhy 

dz 



jci_ sm If, 



Jc± cos if, 



(3) 
(4) 



and one has only two equations for the three functions 
hx{z), hy{z), (p{z). 

In real samples of nonsymmetric shape, adjacent flux 
lines may be slightly rotated relative to each other in 
the critical state. It is this rotation that generates a 
component of the current along the magnetic field. The 
rotation of flux-lines can lead to their mutual cuttingiSii 
Flux line cutting occurs when the component of the cur- 
rent density parallel to the magnetic field, j|| , exceeds 
some longitudinal critical current density ■ In this 
situation a vorte?*^ or a vortex array^, becomes unstable 
with respect to a helical distortion, and the growth of 
this distortion leads to fiux-line cutting. When both j|| 
and jci. are equal to their critical values and jc±, re- 
spectively, Clem's double critical stataStiS occurs in the 
superconductor. In this case one has three conditions for 
the three quantities, and equations (Q) are sufficient to 
describe the double critical state. However, in many real 
situations, jy does not reach , and flux cutting then 
does not occur in the critical state. It is such situations 
that we consider here. In particular, in the above exam- 
ple the projection of the current density on the magnetic 
field (which practically coincides with the z axis) is neg- 
ligible, and flux cutting does not occur. 

We now show how the critical state problem can be 
solved for superconductors of arbitrary shape. Let the 
critical state be known at some moment of time t, i.e., one 
has H — H{r)u{r) inside the superconductor where the 
magnitude of the magnetic field, H, and the unit vector 
v are both known functions of the coordinates r at some 
external magnetic field Ha (t) . The current density j (r) 
in the critical state follows from the Maxwell equation 
j = rot (r)i>'(r)), while the component of the current 
density perpendicular to the magnetic field is given by 
= j — I'ii'i) = jc±n±{r). Here the last equality defines 
the unit vector n^. Let the external field infinitesimally 
(and slowlyii) change by SHa — tlaSt. We now shall 
find the new critical state at the new external magnetic 

field Ua + SUa. 

Under the change of Ha, the critical currents locally 
shift the vortices in the directionpi^ of the Lorentz force 
[j X f]; this shift generates an electric field directed 
along [i' X [j X v]] — i.e., along the vector n_L. 
Thus, we can represent the electric field E(r) in the form 
E = n±e where the scalar function e(r) is the modulus 
of the electric field. Note that in contrast to the Bean 
assumptioufii the electric field generally is not parallel to 



rot(en^) = — /ipH, 
where H = dH/dt, and the equation 



(5) 



(6) 



one can express the change of the magnetic fields and 
currents via one scalar function e(r). This function can 
be found from the condition that in the critical state 
the absolute value of is a given function of B, j± = 
jc±(B), or in the differential form. 



Jx-^=jcx(B)(^^.MoH 



In our case when jc± =const., this condition reads 

Taking into account the definition of j_L and using the 
identities 

Hi> = B.-{u- H)iy, 

J ■ V = {u ■ i-oti/)H, 

we arrive at an equation for e(r), 

■ {rot rot(en^) — {v ■ rotv) rot(en^)} = 0. (7) 

Continuity of the magnetic field on the surface of the 
superconductor, S, yields the boundary condition: 



rot (e(rs)n_L(rs)) = ^J.oi^a + 

[R X rot rot(e(r')ni(r'))] 
47ri?3 



dr\ 



(8) 



where rs is a point on the surface S, R = — r', 
R = |R|, and the integration is carried out over the vol- 
ume of the sample. The right hand side of this bound- 
ary condition expresses ^qH on the surface of the su- 
perconductor (but reaching from outside) with the use 
of the Biot-Savart law. If in the critical state of the 
superconductor there are also boundaries at which the 
direction of the critical currents changes discontinuously 
or which separate regions with j± = jci_ from regions 
with j = Ofi^ the function e(r) has to vanish at these 
boundaries. Otherwise, the electric field en^ would be 
discontinuous there. 

After determining the function e(r), one can find the 
new critical state H(r) -I- SH{r) using the definition 
(5H(r) — H(5t and Eq. We emphasized that the new 
critical state depends only on the previous state H(r) 
and on the change of the external field SHa = Ho^i- 
The dependence on SHa follows from the proportionality 
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of e, H, dj±/dt to Ha, which results from the linearity 
of Eqs. ©-(IHl)- Note that in agreement with the mean- 
ing of the critical state, the new state will be the same 
for different sweep rates of the external magnetic field. 
Ha, since it depends only on the product ilaSt = JHa- 
On the other hand, the electric field e plays an auxiliary 
role in the above description since it is proportional to 
Ha rather than to 6tia- The presented description also 
shows that the critical state generally depends on the his- 
tory of its creation. In other words, it depends not only 
on the final value of the external magnetic field Ha but 
also on the sequence of steps (5Ha that lead to this value. 

The above approach, which is in essence the generaliza- 
tion of the appropriate analyses used for a slabfiSJi^ can 
be summed up as follows: We add to the static equations 
the quasistatic Maxwell equation JSJ. It is knowni 
that for this set of the equations to be solvable, it has to 
be supplemented by some law E(j). We introduce this 
law from well-known physical ideas: At any given j and B 
(determined by the previous critical state), the direction 
of E follows from the formula E = [B x v] where v is the 
vortex velocity caused by the Lorentz force [j x B] . As to 
the magnitude of E, it is found from the condition that 
|j^| = jc_L. In fact, this condition may be interpreted 
as a current -voltage law with |E| = e = at j± < jc± 
and e — > oo at > jc±, which is usually implied in the 
description of the ideal critical stateiii 

Recently^Siii a variational principle was put forward 
to describe the critical state in superconductors. In de- 
riving this principle Badia and Lopez used Eqs. (Q), ® 
and a current-voltage law with \E\ = at j < jc and 
E| — !■ cx) at J > jc- However, the physical idea on the di- 
rection of the electric field was not incorporated in their 
theory. For some situations this leads to contradiction 
with existing concepts, ^'^ In particular, in their so-called 
isotropic model with i7-independent jc, the electric field 
E is parallel to j, and hence a nonzero E along H appears 
even for an infinitesimal longitudinal component of j, i.e., 
flux-cutting occurs without any threshold 

To illustrate the obtained results, we now consider the 
example mentioned above. In the case of the slab, equa- 
tion (JT)) for the electric field e takes the form: 



e"-(^')'e-0. 



(9) 



where (p' = dip/dz and e" = d^e/dz^. For the angle (p 
we obtain from Eqs. (O and (|SJ): 



dip nil II 
MoJc± ^2ep +ep . 
ot 



(10) 



These equations complement Eqs. Q, Q, and now we 
have four equations for four functions. The boundary 
conditions to Eqs. ©, ©, dll) at 2 d/2 are 



{e cos ip)' = 



'Mo 



dhay{t) 

~~dt ' 



hy - Ky{t), (11) 

(esinv3)' = 0, (12) 
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FIG. 1; Profiles of the angle of currents •p(z), magnitude of 
electric field e(z) (top) and magnetic field components hx{z) 
(dashed lines), hy{z) (solid lines) (bottom) in the critical 
states of the slab described by Eqs. ©, EU, l(TT^ - 

inil; /lax = 1.1 and hay = 0.01, 0.04, 0.1, 0.2, 0.35, 0.5, 0.7, 
1, 1.4, 2, 4. We start at hay = with hx{z) = hax — 1 + z, 
hy{z) = 0, and (j}{z) — n/2. Here z is in units d/2, hx and hy 
in units jc±d/2 — Jc/2, and e in units iJ,o(dhay/dt)d/2. 



or equivalently, conditions p2|) . which follow from for- 
mula ijSJ, can be rewritten in the form: 



dhay{t) 

-fiO J- C0Sly9, 



eip 



dhay{t) . 

Mo T — sm ip. 



dt dt 
Taking into account the symmetry of the problem— 



(13) 



it 



is sufficient to solve equations I0J, ll^Jl, (|1U|I in the 
region < z < d/2. Atz = 0, where the direction of the 
currents changes discontinuously, one has the additional 
condition for e, 



e(0) = 0. 



(14) 



Since after switching on hax, the critical currents flow 
in the y direction, we have the following initial state for 
Eq. 



ip{z,t = 0) = 7r/2, 



(15) 



where the moment t = corresponds to the beginning of 
switching on hay. 
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In agreement with the general considerations given 
above, it follows from Eq. (|10|l and e cx dhay/dt that 
the angle is a function of z and hay rather than of z 
and t. In fact, this equation describes (p(z) in the se- 
quence of the critical states developed in the process of 
increasing hay The solution of equations JSJl, Q, 
(|1U|I with conditions 1)11(1 - (|15|l is shown in Fig. 1. In- 
terestingly, when hay increases, the other component of 
the magnetic field, h^, penetrates further into the slab 
[at hay Jc = jc±d, hx{z) almost coincides with hax], 
and the angle if tends to tt. In other words, with in- 
creasing hay the initial critical state for the component 
hx relaxes, while the critical state for hy is developed. 
Note that we should arrive at a different critical state 
with If = tt/2 + &rctAn{hay/hax) if the x and y compo- 
nents of the external field were increased simultaneously 
[hay{t)/hax{t) =const.] from zero to the same values hax, 
hay ~ Jc- Thus, the dependence of the critical state on 
its prehistory is clearly seen even in this simple example. 

Figure 1 also reveals the following two interesting fea- 
tures of the critical state: (a) The visible penetrating 
front of hy reaches the center of the slab when hay is 
still less than Jc/2, the field of full penetration in the 
Bean case, (b) The change of the angle ip{z, hay) has 
diffusive character. This is in stark contrast to the usual 
Bean critical state, in which any change of the current 
direction occurs inside a narrow front. 

Interestingly, equations ®, l(Tn|l are applica- 

ble also to a number of other physical problems if the 
boundary and initial conditions are changed appropri- 
ately. In particular, these equations also describe the 
usual Bean critical state in the slab, corresponding to a 
discontinuous solution ^p{z). Using these equations, one 
can also investigate the low-frequency response of the 



slab to a circularly polarized ac field applied in the plane 
of the sample perpendicularly to the large magnetic field 
Ha^ It is clear from the data of Fig. 1 that this re- 
sponse will differ from the response to a linearly polar- 
ized ac field, for which the analysis based on the usual 
Bean model is applicable. These equations also enable 
one to consider the vortex-shaking effect:^- If the field 
Ha is not uniform in the plane of the slab, and thus a 
sheet current J flows in it, a small ac field applied along 
the current leads to a continuous drift of vortices in the 
direction [J x Hq]. It turns out that Eqs. ||3J), (©, El, 
(fTn|l have a solution that reproduces this result of Ref.l2(]l 
obtained there by a different method using geometrical 
arguments. 

In summary, we have extended the critical-state the- 
ory to the general case when the sample is not sufficiently 
symmetric, or when the external field is not along a sym- 
metry axis or its direction changes in some complex man- 
ner. In such situations the currents in the critical state 
need not be perpendicular to the local magnetic fields, 
and a longitudinal component of the currents with re- 
spect to these fields exists in the superconductor. When 
the magnitude of the longitudinal current density jy does 
not exceed some critical value the critical state 

can be found using our Eqs. (O - (|SJ|. Such a state, in 
general, essentially differs from the usual Bean critical 
state. In other words, the Bean state is only a special 
case of the general critical state. When the component 
j|| reaches jc\\ in the sample (or in part of its volume), 
Clem's double critical state develops thereiiS 
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